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Abstract 



A Lagrangean for the dynamics of an electromagnetic field in a dis- 
persive and dissipative material is constructed (adapting some ideas 
by Bekenstein and Hannay) and an expression for the energy density 
that is positive is obtained from it. The expression contains extra 
(sink) degrees of freedom that represent dissipating modes. In simpli- 
^Sj ■ fied cases the sink modes can be eliminated to yield an energy density 

expression in terms of the electromagnetic fields, the polarization and 
the magnetization only, but which contains parameters associated with 
the sink modes. The method of adding extra modes can be used to 
set up a Lagrangean formalism for dissipative systems in general, such 
that will reinstate time-translation invariance and will yield a unique 



o 
o 

q ■ energy density. 

j^ : 1 PACS: 

Or 

03.50.De; 71.36,+c 
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3 Background 

The problematic nature of the electromagnetic field energy in a dissipative 
material is apparent already at a graduate-teaching level: a frequently used 
textbook notes that in a dispersive medium the energy density lacks a ther- 
modynamical interpretation [1]. (A dispersive material is necessarily dis- 
sipative, since by the Kramers-Kronig relations its constitutive constants, 
that are by definition frequency (to) dependent, must have an imaginary 
part, which represents absorption or energy loss.) The expression for the 
energy density that was derived in [1], was valid only for fields that are 
nearly mono-chromatic. For other basic treatments we quote textbooks as 
[2]- [4], and note that the last reference labels the expression obtained under 
the above conditions as an effective energy density. 

We write the relations between the electric field and the displacement 
using the frequency dependent permittivity e(u>) and between the magnetic 
field and the induction using the permeability h(uj), as 

D(x, uo) = e(uj)E(x, uj), B(x,uj) = /i(uj)H(x,uj). (1) 

Then according to [4] we obtain an effective energy density: 

' d(ue) 



u e f f = Re 
+ Re 



dm 

d(ujfi) 

du 



(<*>o) 



(E(x,t) ■ E(x,t)) 
(H(x,t) ■ H(x,t)) 



The brackets () designate an averaging over a period of — , where loq is the 
carrier frequency. 

Alternative expressions that were subsequently proposed were contro- 
versial. One derivation, that required a significant departure from standard 
electromagnetic theory [5], postulated the independence of the expression of 
particular material constants. This was strongly criticized in [6], where the 
energy-density formula contained explicitly parameters that were present in 
the equation of motion for the field. Another derivation, in terms of consti- 
tutive constants [7], was found to lead to energy densities that are negative 
for a medium with a narrow resonance [8]. 

The quest for a non-controversial energy density can be regarded from a 
more fundamental angle. This quantity is expected to be part of a conserva- 
tion equation, which (by a well-known theorem due to Noether) is tied to the 
time invariance of the Lagrangean. As soon as the invariance is lost (which 



is the case for dissipative systems), the quantity to be conserved is unde- 
fined. Thus the need to find an appropriate energy density is expected to 
arise in a wider context, too. (In a dramatic account, in which the proposer 
of a time-varying light- velocity describes his Iliad to get his idea published, 
one also finds the problem of the proper Lagrangean formulation to play a 
key role [9].) 

The issue of the electromagnetic-field energy density in a dispersive and 
dissipative medium has recently resurfaced in the contexts of the sublu- 
minality of light-propagation and of left handed materials ([10], [11]). The 
former work utilizes some analytical properties of the constitutive relations, 
while the departure point in the latter are the equations of motions for the 
macroscopic polarizabilities. As an application, the energy density for a left 
handed medium has been calculated in [11]. The expressions in the two 
papers differ. 

Our approach is based on the recognition that in standard field theories 
(classical or quantal) the energy momentum tensor satisfies a conservation 
equation having the form of a "div" -equation involving all degrees of free- 
dom. Thus, it seems, that one should be able to write out an energy-density 
in an unambiguous fashion starting with a Lagrangean. (The idea of a La- 
grangean formulation was raised earlier [12], but only in a programmatic 
manner and by making an assumption that the detailed treatment worked 
out in this paper does not justify.) 

As is well known, the energy density Tg is a component of the energy 
momentum tensor T- , which is uniquely derivable from the Lagrangean 
density C(uk,Ukj), this being a function of the field variables u^ and of 
their derivatives u^j- (The indexes take the value for the time component 
and the numbers (1,2,3) for the remaining, space components. Summation 
is implied for repeated indexes and a symbol after the comma represents 
a derivative. The dot will also be used for a time derivative.) The formal 
definition is given by 

T* = ^u^-a) (3) 

Two circumstances appear to prevent one from deriving the energy- 
density from the conservation equation in a dissipative material. First, a 
Lagrangean has not been formulated for the equations of motion (but this is 
done below). Secondly, the energy sinks involved in the dissipative mecha- 
nisms have not been given a dynamic representation, but only a phenomeno- 
logical one (see equation (1) ), namely, through the appearances of complex 



permittivity and permeability in the constitutive relations or, equivalently, 
through the presence of time-reversal non-invariant terms in the equations 
of motion. 

The first issue was recently resolved in a short note [13], where it was 
shown how to formulate the Lagrangean (actually, the Hamiltonian) for a 
dissipative case. The second point was (indirectly) addressed in [14], fol- 
lowing a previous publication in [15]. These works formulated the dynamics 
of a (conjectural) time dependent fine-structure constant by including in 
the Universe an additional degree of freedom. We follow up both these 
approaches, with changes required by the different context. 

4 A Lagrangean for static polarization and mag- 
netization 

Before approaching the general problem we will first present a Lagrangean 
density for the electromagnetic field in a material which has a static polar- 
ization and magnetization. The following expression connects the electric 
displacement field D with the electric field E and polarization P: 

D = eE + P (4) 

(e is not necessarily the vacuum permittivity eo). Similarly the magnetic 
field H is connected with the magnetic induction field B and magnetization 
M by: 

H=-B-M (5) 

I- 1 

(again ji is not necessarily the vacuum inverse permeability jjlq). The above 
fields satisfy both the homogeneous Maxwell's equations: 

curlE + B = 0, divB = (6) 

and the inhomogeneous equations: 

curlH - D = J, divD = p (7) 

In the following we assume that both free charges p and currents J are 
zero. The above equations can not be obtained from a Lagrangean density 
expressed in terms of those fields. However, this situation can be amended 



by representing the fields in terms of vector A and scalar $ potentials, as 
follows: 

E = -V$ -A, B = curlA (8) 

Using these definitions we see that the homogenous equations (6) are satis- 
fied automatically. The inhomogeneous equations (7) can be obtained from 
the functional derivative of the Lagrangean: 

£ = £eb + £pe + C'MH (9) 

Ceb = \[tE 2 -^] (10) 

C PE = PE (11) 

C M H = liM-H (12) 

in which the reader should think of all field quantities as given in terms 
of the potentials of equation (8) . Using equation (3) one can obtain the 
energy density: 

fir 

T ° = ?t.A-£. (13) 

dl 

which yields for the Lagrangean of equation (12) after some manipulations 
the result: 

T ° = i [eE 2 + fiH 2 + fiM 2 ] (14) 

The last term juM 2 can be usually discarded for static magnetization since it 
represents nothing but a time-independent constant. However, we include it 
for future sections in which the polarization and magnetization will become 
dynamic degrees of freedom. 



5 Equations of motion for polarization and mag- 
netization 

In cases that the polarization P is induced by an electric field E, and the 
magnetization M is induced by a magnetic field H, one needs a set of equa- 
tions to describe these processes. The equations are given, e.g., in [11] as: 

P + T e P + u 2 P = eo^E (15) 

M + T h M + uj\M = FujIH (16) 



with all symbols and notation as denned in [11]. In addition, one has the 
Maxwell-equations given in the previous sections and derived from the La- 
grangean of equation (12) using the contravariant four- vectors A a = (A, 3>). 
Turning now to dissipative cases, following [13] -[15] we introduce a set 
of new fields, here named canonical fields, which are formally distinguished 
from the previous, physical fields by writing them in lower case symbols. 
These fields (assumed to be real) are made to be part of a Hamiltonian 
or Lagrangean formulation, in contrast to the preceding physical quantities 
(written in upper case symbols), with which this cannot be done [13]. Ex- 
plicitly, we shall work with the scaled polarization p a and magnetization 
m a ; the electromagnetic fields e a , h a , the (contravariant) vector-potentials 
a a , and (again following [13]) also the fields r a ,n a conjugate to p a and to 
m a (as defined below). The choice of the fields (and of the various constant, 
scaling factors) is guided by the requirement that we recapture the form of 
the equations of motion and the energy density currently widely employed 
in solid state optics in the appropriate limiting case of no dissipation. (This 
requirement is at times detrimental to the formal symmetry between the 
modes linked with the polarization and the magnetization, respectively.) 
The physical meaning of the canonical fields will be made clear later by 
comparing their equations of motion with those of the physical fields. (For 
a similar procedure see [13].) In addition, we shall introduce two new scalar 
fields: first \I/, representing a degree of freedom associated with the dissipa- 
tion of the polarization (analogous to the degree of freedom denoted with 
the same symbol in [15]) and, secondly, the "sink-field" $ for the dissipation 
in the magnetic mode. 

6 The dissipative Lagrangean 

This has the following parts: 

L, = L, e \) + L,p e + J-"mh ~r L-pr t L-rnn ~r J~,\j/ + £,<j> + Lsp r —\$i + L. mn —Q \l I ) 

This Lagrangean contains first the Lagrangean of equation (12) , there writ- 
ten in terms of the electromagnetic fields, and now re-expressed in the new 
lower-case variables as C e b + C pe + £ m h- In addition, the above expres- 
sion contains the Lagrangeans in the polarization modes, the magnetization 
modes, the electromagnetic fields, the ^ and <5 sink fields and then the 
interaction-Lagrangeans between polarization and electric field, magnetiza- 
tion and magnetic field and, finally, the polarization and magnetization with 



their respective sink fields. Writing these out fully: 
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(26) 



having chosen the polarization and magnetization vectors to be along the 
x (or 1) axis. Physically, the new degrees of freedom \1/ and $ can be 
associated with some sort of relaxation mechanism for the polarization and 
the magnetization, respectively. In the above, we have chosen for Cq and 
£$ what appear to us the simplest type of expressions that can represent 
unspecified (and, so far, arbitrary) degrees of freedom. The symbols k, p, A, v 
stand for constant, non-negative coefficients, whose values depend on the 
nature of the relaxation mechanisms. 

Equating to zero the variation of £ with respect to each degree of freedom 
leads to the equations of motion, in accordance with the Euler-Lagrange 
equations. For the variables r and n the equations obtained are: 



r 
n 



P- 
fn 



- m& 



(27) 
(28) 



which can be inserted into the equations for p and m. These take now the 
following form: 



p 



m 



(§ + ijr 2 _ w 2^ 

($ + 6 2 - uol)p 



eo^ei 
FojQhi 



(29) 
(30) 
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The electromagnetic equations of motion are identical to Maxwell's equa- 
tions in the new (small letter) variables. Finally, for the sink coordinates 
one has 

k Ateo^p at 

■■ v T u d(mn) , . 

#+ A* = ^"W i (32) 

From equation (3) we obtain the expression for the energy density Tq , 
namely, 

tf = E5|*-* < 33 > 

where the sum is over all degrees of freedom. We separate the electromag- 
netic, the polarization and the magnetization parts 

T o = ( T o )e,h + (T )p,v + (T ) m ,* (34) 

where the first term has the well known form of the electromagnetic energy 
density given in equation (14) and repeated here for completeness: 

(T °k, = f e 2 + f ft 2 + f rJ (35) 

For the others we obtain 

(T °W = ^2^-^) 2 +^V] + f* 2 + f* 2 (36) 

X (p 2 + c r V) + ^ 2 + gvi/ 2 + -l^p 2 vi/ 2 



2e u;2 - ' - " ' ' 2 " ' 2 " ' 2e w* 
1 9p 2 



2eo^p dt 



* (37) 



and 



2 -|^[ ( m-m<i>) 2 + ^ 2 m 2 ] + ^ + ^ 



WW - ^[(™-™ $ ) +^™1+o* +o* ( 38 ) 



MO , . 2 i 2 2\ , ^i.2 i ^^2 i ^0 2A2 

2Fw 2 v ° ' 2 2 2Fw 2 

Mo dm 2 ■ 

m - $ 39 

2Fw 2 9t v 7 



The first-written form (=sum of squares with non-negative coefficients) of 
either quantity guarantees that each part of the energy density is positive 
(non-negative). The usual expressions for the energy density, e.g. in [4] or 
[11], differ from the above by the presence of the terms in <J> and \l/ and 
their time derivatives. In the following section it is our purpose to eliminate 
these variables by making use of the equations of motion, equation (31) and 
equation (32) . 

7 A solvable case 

We take a simplified case when the coefficients p and v are both zero. Then 
the equation of motion for the polarization sink-mode equation (31) can be 
integrated. (A similar procedure applies to the magnetization sink- variable.) 
We assume the following initial conditions for *&(t) 

*(0) = 0, tf (0) = T e /2 (40) 

The reason for these choices is that with them for short times, t << — the 
differential equation (15) is regained. (This will be shown presently. Cf. 
also [13].) Then from equation (31) 

# = ^_-gv( t) - l 9 -&m ( 4i) 

2neo(jj 2 dt 2 Keo^p dt 

Integrating once 

1 dp 2 (t) p 2 (t)ty „ , s 

ZneoLj* ot neocu- 

(l + i%* = > Mfi + C (43) 

K€qC0~ IKEqUJ^ Ot 

with the (first) integration constant given by 

r e(1 | p 2 (Q) ) p(o)p(o ) 

2 neoLO 2 hequj 2 

Integrating once more and arranging for satisfaction of the first initial con- 
dition in equation (40) , we finally obtain: 



leading to 



C= i£ ( l + ^_)_^^ (44) 



*(t) = -In 



, (l + «, 



neou)* > 



2"' L (1+ Pi(o^ 



r e P 2 (o) p(o)p(o) 1 /•* df 

+ M 1 + -Z— 2) " ~ ..2 1 /. : ^777V ( 45 ) 



When this expression and equation (42) are substituted into equation (37) , 
one obtains after considerable simplification the following expression for the 
energy density arising from the time varying "canonical" polarization p(t): 

(^ 1 ( p 2 (t) + Ke^ 2 p C 2 22 s 

When we assume that T e is small (precisely, T e ^fH « 1), and so is p(0) 
(this will be confirmed in the next section), then C 2 will be a second order 
correction which may be neglected: 

^-JS?(lfsf +, *"w) (47) 

A similar expression is obtained for the part of the energy density involving 
the "canonical" magnetization m(t), the variables n(t) and <J>(£) having been 
eliminated through their equations of motion, 

(T ° )m '* = 2KI ( 1 T Wit) + w ° m (t)) (48) 

U "•" AFc^ 

The above expressions, equation (47) and equation (48) , are quite similar 
to those in equation (11) of Ruppin [11], except that they are written in the 
canonical (small letter) variables, rather than in the physical variables (for 
relations between these, see immediately below), and that they contain time 
dependent denominators. 

The main results of this work, equation (46) and its analogue for the 
magnetization energy density, are exact and contain nonperturbative cor- 
rections to the energy density, due to the presence of the sink degrees of 
freedom. While exact, they are model dependent in the sense that sinks 
represented by different Lagrangeans would lead to different energy densi- 
ties. This is clear, due to the presence in the energy densities of the param- 
eters k and A that were introduced in the Lagrangean in equation (23) and 
equation (24) . It is of interest to note that the non-dissipative limit is not 
regained when T e , T^ — ► 0, but only when also k, A — ► 00. 
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8 The physical fields 

To recapture the basic equations of motion equation (15) and equation (16) 
for the physical polarization variable, we proceed as follows: We postulate 

p(t) = e*(*)P(t) (49) 

m(t) = e*®M(t) (50) 

Then, from equation (40) , for short times < t << p-, 

* « *(0) + *(0)£ = r e i/2 (51) 

This turns (the vector form of) equation (29) into the following: 



r e t 



P + T e P + co 2 r P = euj 2 p e- — e (52) 

and likewise for the magnetization variables. Recalling equation (15) and 
equation (16) , we can thus extrapolate to later times so as to identify 

e(t) = e*®E(t) (53) 

rn(t) = e^H{t) (54) 

Thus, all the "canonical" variables are the fields in which the decay of the 
physical field variables has been reinstated. On the other hand, the decay is 
itself dependent on the fields. (Cf. [12].) Furthermore, Maxwell's equations 
for the physical fields are also modified, just as in [14] and [15] (eq. (14) and 
cq. (12), respectively). 

9 Conclusion 

Using standard Lagrangean formalism, we have obtained a unique energy 
density for a dissipative medium (capable also of sustaining an electric po- 
larization and/or a magnetic one). It has been found necessary to introduce 
two additional degrees of freedom ("sink-modes"), associated with decay 
mechanisms in the electrical and magnetic modes. The derived contribu- 
tions to the energy density, shown in equation (46) - equation (48) , are all 
positive and are similar in their form to the corresponding results in [11]. 
They contain particular physical parameters, as anticipated in [6]. The 
energy density is part of a conservation equation involving also an energy- 
current (momentum) density (the Poynting vector) formally given by the 
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tensor-components T$(a = 1,2,3) shown in equation (3) . These properties 
appear to be true in general, namely, for a large variety of systems that are 
dissipative. 

A solvable model has been worked out in this paper. Both in this case and 
in a general one, the energy density has corrections not appearing in other 
approaches not based on a Lagrangean. On the other hand, in a spatially 
homogeneous infinite medium (such as treated in this paper) the momentum 
density is still the Poynting vector. This can be seen from equation (3) , 
since the second term in that expression is absent (for a differing from 
(3) and there are no space derivatives in the Lagrangean, other than in 
the electromagnetic part. (As a consequence, the Lagrangean is not fully 
Lorentz- invariant . ) 

For a medium possessing space varying properties (including, under cer- 
tain circumstances, a finiteness of size) the procedure outlined in this paper 
can be extended to space-coordinates. This means the introduction of new, 
space-dependent "sink" degrees of motion (e.g., at the boundaries). As a 
consequence, additional terms would appear also in the momentum density. 
Apart from an interesting extension of the Poynting vector concept, the 
suggested approach can also have the practical use of solving electromag- 
netic problems involving dissipation, such as laser gain and loss calculations, 
microwave losses in wave guides and cavities and transformer performance. 
The approach can be further extended to other fields of macroscopic physics 
such as viscous fluid dynamics. Work in these directions is under progress. 
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